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The interaction of optical vortices (or phase singularities, screw dislocations) with ordinary matter is treated 
with simple approach. Using total internal reflection phenomenon and superposition of four plane waves 
incident on a material with refractive index lower than the original propagating medium, we are able to show 
birth and annihilation of optical vortices in the evanescent field with curved topological features. Until now, 
this phenomenon has been explored only in free space propagation. By a suitable tuning process, involving 
the incident angles and the amplitudes of the incident plane waves, it is possible to create unusual topological 
features of optical vortices in the vicinity of the material. We believe that this work can open new aspects of 
curved optical vortex manipulation in near-field optics. © 2011 Optical Society of America 
OCIS codes: 260.0260, 260.6970, 050.0050, 050.4865. 



Optical vortices (or phase singularities, screw disloca- 
tions) are now well established in the scientific commu- 
nity as basic entities of the electromagnetic field. The 
first work on optical singularities was carried out by Nye 
and Berry [1] in the 70's. Almost twelve years later, sin- 
gularities were recognized closely related to the orbital 
angular momentum of a beam of light [2] and now we 
can say that classical optics is able to describe very well 
orbital and spin [3] optical angular momentum contribu- 
tions in a beam of light. 

If the phase around a point in the field increases (or 
decreases) by an integer multiple of 2tt, also called topo- 
logical charge, the phase at the point is undefined. Con- 
sequently the field amplitude drops to zero and we have a 
singularity. The dynamics of topological charges in prop- 
agating fields was studied by many authors [4-11]. A 
useful visualization of the optical vortex is to assign a 
"path" to the point where the phase is undefined, as a 
function of position. In this way we can construct graphs 
which show the dark lines representing the points of sin- 
gularities in the field. For example, the well-know Bessel 
beam [12] has an on-axis singularity because the phase 
dependence of the field is given by (f>{6, z) = f3z + raO 
where (3 is the wave vector in the direction of propaga- 
tion, 9 is the azimuthal angle in cylindrical coordinates 
and m is the topological charge of the beam. If we plot 
the function (j){9, z), we observe that a singularity exists 
(for m 7^ 0) in the center of the beam [p = 0) and follows 
a straight line propagation parallel to the z-axis. Also, 
the surface where the phase is constant (i.e. the wave- 
front of the beam) is helical. Such helical wavefront give 
a non-zero contribution to the transversal component of 
linear momentum. This can be used to move particles 
around the singularities as opposed to the spin angular 
momentum which rotates the particles around their own 
axis [10]. 

The problem of producing optical beams with vortices 
paths that are not straight lines, arises naturally and 



was answered theoretically by M. V. Berry and M. R. 
Dennis [13] by making interference of Bessel beams and 
using a perturbative analysis. The theory of Berry and 
Dennis was verified experimentally by Leach et. al. [14]. 
In particular, O'Holleran et. al. [15] showed that we can 
have unusual topological features (loops, births and an- 
nihilations) arising with the coherent superposition of 
three, four and five plane waves. Their work inspired the 
present paper because a rich variety of singular paths can 
exist, using only these basic entities of the electromag- 
netic theory, i.e., plane waves. Non-straight-line propaga- 
tion of optical vortices is also found in free-propagating 
Airy beams, which results in parabolic trajectories, as 
demonstrated by Mazilu et. al. and Dai. et. al. [16,17] 

In this work, we used the total internal reflection phe- 
nomenon to create evanescent fields with curved topo- 
logical singularities. The papers cited previously lead us 
to consider questions like if is it possible to have loops, 
knots, birth and annihilation of optical vortices in the 
evanescent field. Also, the question of what happens to 
the optical vortices when they interact with matter is 
studied. Here we show that it is possible to visualize, 
numerically, birth and annihilation of singularities. The 
question involving loops and knots remains open to nu- 
merical and experimental investigation. 

In order to show the relevant aspects of the effects 
in near-field regime, we take a planar surface separat- 
ing two media with different index of refraction. Plane 
waves (four, in this paper) with wave vectors k trav- 
elling originally in the medium with index of refrac- 
tion rii making an angle 6j with the normal, strikes a 
medium with index of refraction n2 giving rise to re- 
fracted waves with wave vectors k" and reflected waves 
with wave vectors k' . If we take ni > n2, there exists 
an angle 6c ~ sin^^ (77,1/712), namely, the critical an- 
gle, such that if the incident plane waves satisfy the 
condition 6j > 9c the cosine of the refracted angles 
becomes purely imaginary, consequently, the waves re- 
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main bounded within the z direction with dependence 
exp{—z/d), where d^^ = k'{sin^9j — sin^Oc)^/'^ is the 
characteristic length scale. 

The evanescent field in the region 2; > can be calcu- 
lated by [18] 



Eev{x,y,z > 0) 



4 

E 



2niEjo cos 6j 



j^l J ~<~ V"2 ~ "1 """^ 



rii cos 

X exp[ik"{x cos (pj sin 9^)] 

X exp[ik" {y sin 9 j sincpj + 0cos6'^)] 



(1) 



where 9j and 9j arc the refracted and incident angle of 
the j-th incident wave, respectively, Ejo is the amplitude 
of the j-th incident plane wave and 4'j is the azimuthal 
angle in spherical coordinates. Using Eq. (1) and varying 
twelve parameters, the incident angles and the incident 
amplitudes, we can calculate the evanescent field above 
the plane for a variety of situations. 

Let us take four plane waves propagating from medium 
with refraction index ni = 1.5 to the medium with re- 
fraction index n2 = 1.0. We assume that the polar an- 
gles in spherical coordinates are given by 9i = 1.29c, 
6*2 = 1.50c J ^3 ~ 1-66'c; 94 = 1.79c, the azimuthal an- 
gles are given by 01 = 0, 02 = n/2, ^3 ~ n, (j)^ — 3tt/2 
and all four incident plane waves have unit amplitude. In 
[15] it was shown that if we superimpose four plane waves 
with equal magnitudes in free space, the singularities are 
straight lines. So, this is what we get for z < before the 
waves strike the surface. But when we look at the field 
right above the surface, we observe that there are line 
vortices annihilated, as shown in Fig. 1. The presence 
of the boundary, which separates two different media, 
breaks the symmetry of straight line vortex propagation. 
To explain this result, we should look at the amplitudes 
of the individual refracted waves. They all depend on z, 
so wc have plotted in Fig. 2(a) the dependence of ai + 04 
and 02 -l- 03 {uj is the amplitude of the j-th refracted 
plane wave) as a function of z. Wc clearly see that for 
any value of z above the surface, the following inequal- 
ity is satisfied: ai 4- 04 > 02-1-03. As demonstrated in 
[15] this condition implies that loops of vortices must be 
present. The Fresnel's coefficients play the principal role 
in this discussion because they are the amplitudes of the 
refracted waves necessary to account for the phasor de- 
scription (see [15]). Also, we observe that they appear 
in pairs, as expected for topological charge conservation 
[19]. Finally, part (b) of Fig. 1 shows the path of the vor- 
tices marked with a red and blue circle as a function of x, 
y and z. Clearly we see the annihilation effect of vortices 
appearing solely due to the spatial symmetry breaking. 

Creation of pairs of singularities is also possible. In 
Fig. 3(a) we show the phase profile of the evanescent 
waves but with different parameters given by 9i = 92 = 
9^ = 1.29c, 04 ~ l.T9c, the azimuthal angles are the same 
as in Fig. 1. and the amplitude of the fourth incident 
plane wave is taken to be 5.5. With these parameters the 




Fig. 1. (Color online) Annihilation of optical vortices in 
the evanescent field. Parameters used in the simulation: 
9i = 1.26ic, 6*2 = 1.56*^, 6*3 = 1.69c,94: = 1.761^, 0i = 0, 
4>2 = 7r/2, 03 ~ n, (j>4 ~ 37r/2, ni = 1.5 and 712 = 1.0. 
(a) Phase profile of the evanescent field. The blue and 
red lines indicates singularities with opposed topological 
charges, (b) The vortex path of the singularities indi- 
cated in part (a) as a function of z. 



amplitudes of the refracted waves are plotted in part (b) 
of Fig. 2, so we still have the condition ai -|- 04 > 02 -I- 03 
satisfied for all z (we were labeling the amplitudes in in- 
creasing order, but in this case there is no difference). 
But now we can demonstrate the birth effect of singu- 
larities with opposed signs, as shown in Fig. 3(b). In the 
picture of Talbot cell, we are simply moving the cell to 
contemplate another region which shows the effect. A 
detailed analysis for the conditions under which we have 
creation/annihilation of optical vortices is still open to 
discussion. The difficulty arises because for each aj we 
have three free parameters {9j,(t>j and i?jo) that, when 
placed in 04-1-04 > 02 -I- 03, generates an inequality with 
twelve free parameters (fifteen if we consider the relative 
initial phases). We are by no means in a position to dis- 
cuss the general case in the present letter, leaving this 
interesting analysis for future work. 

Although we have showed that it is possible under 
suitable tuning process to create and annihilate optical 
vortices in the evanescent field, the symmetry also breaks 
even in the condition 9j < 9c, i.e. for propagating fields 
beyond the boundary (or for ni < n2, in which there is 
no condition for a bounded wave). We choose to adopt 
the evanescent-view for its implications in near-field op- 
tics, still, to our knowledge, unverified in experimental 
conditions. Some questions still need to be answered and 
experiments need to be performed. Also, different mate- 
rials should lead to different and unexpected curved vor- 
tex topological behavior. We do not face these problems 
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Fig. 2. (Color online) (a) Amplitudes of ai + 04 and 
0.2 + 03 in function of z for the annihilation of optical 
vortices case, (b) Amplitudes of a\ + and a-i + 03 in 
function of z for the birth of optical vortices case. 




Fig. 3. (Color online) Creation of optical vortices in 
the evanescent field. Parameters used in the simulation: 

01 = 1.20,, 02 = 1.20c, 03 = 1.20c, 04 = 1.70c, (/-i = 0, 

02 = 7r/2, 03 = TT, 04 = 37r/2, ni = 1.5, n2 = 1.0 and 
the amplitude of the fourth plane wave is 5.5. (a) Phase 
profile of the evanescent field. The blue and red lines 
indicates singularities with opposed topological charges, 
(b) The vortex path of the singularities indicated in part 
(a) as a function of z. 



in this paper, retaining our results only to show that it is 
possible to create new degrees of freedom in the optical 
tweezers area using optical singularities. 

In conclusion we showed that evanescent fields can 
contemplate birth and annihilation of optical singulari- 
ties by simple tuning process. From a physics standpoint, 
the symmetry of straight line propagation is broken when 
they interact with matter, and we showed that the pha- 
sor argument is sufficient to characterize the dynamic be- 
havior of the singularities. As reported by some authors 
[20-22], optical trapping with singularities is a promis- 
ing branch of physics, and this paper on near-field optics 
with curved topological features can open new aspects 
of experimental investigation. 

The authors thanks J. M. Hickmann and S. B. Caval- 
canti for useful discussions. 

References 

1. J. F. Nye and M. V. Berry, Proc. R. Soc. Lond. A 336, 
165 (1974). 

2. L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, and J. 
R Woerdman, Rhys. Rev. A 45, 8185 (1992). 

3. R. A. Beth, Phys. Rev. 50, 115 (1936). 

4. I. V. Basistiy, M. S. Soskin, and M.V. Vasnetsov, Opt. 
Comm. 119, 604 (1995). 

5. J. Courtial, R. Zambrini, M. R. Dennis, and M. Vas- 
netsov, Opt. Exp. 14, 938 (2006). 



6. Y. Zhao, D. Shapiro, D. Mcgloin, D. T. Chin, and S. 
Marchesini, Opt. Exp. 17, 23316 (2009). 

7. J. Leach, E. Yao, and M. J. Padgett, New J. Phys. 6, 71 
(2004). 

8. J. Courtial, and K. OHoIleran, Eur. Phys. J. Special 
Topics. 145, 35 (2007). 

9. N. R. Heckenberg, R. McDuff, C. P Smith, and A. G. 
White, Opt. Lett. 17, 221 (1992). 

10. A. T. O'Neil, I. MacVicar, L. Allen and M. J. Padgett, 
Phys. Rev. Lett. 88, 053601 (2002). 

11. D. G. Crier, and Y. Roichman, Appl. Opt. 45, 880 
(2006). 

12. J. Arlt, and K. Dholakia, Opt. Comm. 177, 297 (2000). 

13. M. V. Berry, and M. R. Dennis, Proc. R. Soc. Lond. A 
457, 2251 (2001). 

14. J. Leach, M. R. Dennis, J. Courtial, and M. J. Padgett, 
Nature 432, 165 (2004). 

15. K. O'HoUeran, M. J. Padgett, and M. R. Dennis, Opt. 
Exp. 14, 3039 (2006). 

16. M. Mazilu, J. BaumgartI, T. Cimr, and K. Dholakia, 
Proc. SPIE 7430, 74300C (2009). 

17. H. T. Dai, Y. J. Liu, D. Luo, and X. W. Sun, Opt. Lett. 
35, 4075 (2010) 

18. J. D. Jackson, Classical electrodynamics (John Wiley 
and Sons, Inc, 1999). 

19. M. S. Soskin, V. N. Gorshkov, M. V. Vasnetsov, J. T. 
Malos, and N. R. Heckenberg, Phys. Rev. A 56, 4064 
(1997). 

20. K. T. Cahagan, and G. A. Swartzlander, Opt. Lett. 21, 



3 



827 (1996). 

21. J. E. Curtis, B. A. Koss, and D. G. Grier, Opt. Comm. 
207, 169 (2002). 

22. K. T. Gahagan, and G. A. Swartzlander, J. Opt. Soc. 
Am. B 15, 524 (1998). 



4 



References 

1. J. F. Nye and M. V. Berry, "Dislocations in wave trains" , 
Proc. R. Soc. Lond. A 336, 165 (1974). 

2. L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, and J. 
P. Woerdman, "Orbital angular momentum of light and 
the transformation of Laguerre-Gaussian laser modes", 
Phys. Rev. A 45, 8185 (1992). 

3. R. A. Beth, "Mechanical detection and measurement of 
the angular momentum of light", Phys. Rev. 50, 115 
(1936). 

4. I. V. Basistiy, M. S. Soskin, and M.V. Vasnetsov, "Op- 
tical wavefront dislocations and their properties". Opt. 
Comm. 119, 604 (1995). 

5. J. Courtial, R. Zambrini, M. R. Dennis, and M. Vas- 
netsov, "Angular momentum of optical vortex arrays". 
Opt. Exp. 14, 938 (2006). 

6. Y. Zhao, D. Shapiro, D. Mcgioin, D. T. Chiu, and S. 
Marchesini, "Direct observation of the transfer of orbital 
angular momentum to metal particles from a focused cir- 
cularly polarized Gaussian beam", Opt. Exp. 17, 23316 
(2009). 

7. J. Leach, E. Yao, and M. J. Padgett, "Observation of the 
vortex structure of a non-integer vortex beam". New J. 
Phys. 6, 71 (2004). 

8. J. Courtial, and K. OHoUeran, "Experiments with 
twisted light", Eur. Phys. J. Special Topics. 145, 35 
(2007). 

9. N. R. Heckenberg, R. McDuff, C. P Smith, and A. 
G. White, "Generation of optical phase singularities 
by computer-generated holograms", Opt. Lett. 17, 221 
(1992). 

10. A. T. O'Neil, I. MacVicar, L. Allen and M. J. Padgett, 
"Intrinsic and extrinsic nature of the orbital angular mo- 
mentum of a light beam", Phys. Rev. Lett. 88, 053601 
(2002). 

11. D. G. Crier, and Y. Roichman, "Holographic optical 
trapping", Appl. Opt. 45, 880 (2006). 

12. J. Arlt, and K. Dholakia, "Generation of high-order 
Bessel beams by use of an axicon" Opt. Comm. 177, 
297 (2000). 

13. M. V. Berry, and M. R. Dennis, "Knotted and linked 
phase singularities in monochromatic waves", Proc. R. 
Soc. Lond. A 457, 2251 (2001). 

14. J. Leach, M. R. Dennis, J. Courtial, and M. J. Padgett, 
"Knotted threads of darkness", Nature 432, 165 (2004). 

15. K. O'HoUeran, M. J. Padgett, and M. R. Dennis, "Topol- 
ogy of optical vortex lines formed by the interference of 
three, four, and five plane waves". Opt. Exp. 14, 3039 
(2006). 

16. M. Mazilu, J. Baumgartl, T. Cimr, and K. Dholakia, 
"Accelerating vortices in Airy beams " . Proc. SPIE 7430, 
74300C (2009). 

17. H. T. Dai, Y. J. Liu, D. Luo, and X. W. Sun, "Propaga- 
tion dynamics of an optical vortex imposed on an Airy 
beam". Opt. Lett. 35, 4075 (2010) 

18. J. D. Jackson, Classical electrodynamics (John Wiley 
and Sons, Inc, 1999). 

19. M. S. Soskin, V. N. Gorshkov, M. V. Vasnetsov, J. T. 
Malos, and N. R. Heckenberg, "Topological charge and 
angular momentum of light beams carrying optical vor- 
tices", Phys. Rev. A 56, 4064 (1997). 



20. K. T. Cahagan, and G. A. Swartzlander, "Optical vortex 
trapping of particles" , Opt. Lett. 21, 827 (1996). 

21. J. E. Curtis, B. A. Koss, and D. G. Crier, "Dynamic 
holographic optical tweezers ", Opt. Comm. 207, 169 
(2002). 

22. K. T. Cahagan, and G. A. Swartzlander, "Trapping of 
low-index microparticles in an optical vortex", J. Opt. 
Soc. Am. B 15, 524 (1998). 



5 



